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Abstract 

A  theory  for  fast  modulation  effects  in  cesium  frequency  standards  has  been  devel¬ 
oped  and  expressions  derived  for  frequency  offsets  produced  by  Rabi  and  cavity  pulling 
with  square-wave  phase  modulation.  These  arc  compared  with  equivalent  expressions 
for  sine-wave  phase  modulation.  Ratios  of  the  amplitudes  of  the  frequency  offsets  in 
the  two  modulation  schemes  are  calculated  for  a  mono-velocity  beam  and  shown  to  be 
in  agreement  with  measurements  performed  on  a  commercial  standard.  This  analysis 
allows  a  more  complete  evaluation  and  comparison  of  different  modulation  schemes 
for  cesium  atomic  beam  standards. 

INTRODUCTION 

A  number  of  different  modulation  schemes  are  used  to  implement  the  servo  loop  in  cesium  atomic 
clocks.  They  include  sine-wave  phase  (frequency)  modulation,  square-wave  phase  modulation  and 
square-wave  frequency  modulation.  The  choice  of  modulation  scheme  is  based  on  a  number  of  consid¬ 
erations,  e.g.,  ease  of  implementation,  optimization  of  cesium  flux  usage  and  Allan  variance,  minimized 
frequency  offsets.  Analysis  of  servo  loop  gain  has  been  performed  for  all  these  schemes  and  expres¬ 
sions  exist  to  model  the  gain  dependence  upon  loop  parameters.  However,  the  analysis  of  the  effects  of 
.modulation  scheme  on  accuracy  and  long  term  stability  are  less  complete.  Sine-wave  and  square-wave 
frequency  modulation  effects  have  been  considered!1!  within  the  context  of  slow  modulation!2!,  hut  fast 
modulation  transients  have  not  been  modeled  because  of  the  lack  of  an  atomic  description  of  those 
atoms  that  experience  the  step  change  in  microwave  phase  or  frequency.  A  theory  for  square-wave 
phase  modulation  is  outlined  in  this  paper. 

In  the  analysis  of  servo  loop  gain,  it  can  be  assumed  that  the  cesium  atoms  experience  a  constant  mi¬ 
crowave  frequency  and  the  effects  of  phase  modulation  only  perturb  the  relative  phase  of  the  precessing 
atoms  and  the  microwave  field  in  the  second  arm  of  the  Ramsey  cavity!3!.  ln  particular,  the  transient 
effects  in  square-wave  phase  modulation  can  be  safely  omitted,  and  the  contribution  neglected  from 
those  atoms  which  experience  the  step  change  in  phase  during  their  resonant  interaction  with  the 
microwave  field.  The  error  incurred  in  the  calculation  of  servo  gain  slope  is  only  of  order  where 
t  is  the  interaction  time  and  is  the  modulation  frequency.  Similarly,  by  neglecting  the  frequency 
variation  of  the  microwave  field,  small  errors  of  order  Awr  are  incurred  in  the  frequency  modulation 
calculations.  ±Acu  represents  the  depth  of  frequency  modulation. 

Spurious  servo  signals  are  introduced  by  neighboring  transitions  (Rabi  pulling)  and  cavity  mistuning 
(cavity  pulling)  and  affect  the  performance  of  cesium  frequency  standards.  Differences  in  the  velocity 


401 


Report  Documentation  Page 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 


1.  REPORT  DATE 

NOV  1989 


2.  REPORT  TYPE 


3.  DATES  COVERED 

00-00-1989  to  00-00-1989 


4.  TITLE  AND  SUBTITLE 

Modeling  Fast  Modulation  Effects  In  Cesium  Atomic  Clocks 


6.  AUTHOR(S) 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

The  Aerospace  Corporation, Electronics  Research  Laboratory,?. O.  Box 
92957, Los  Angeles, CA, 90009 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


10.  SPONSOR/MONITOR'S  ACRONYM(S) 


11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

See  also  ADA224769.  Proceedings  of  the  Twenty-first  Annual  Precise  Time  and  Time  Interval  (PTTI) 
Applications  and  Planning  Meeting,  Redondo  Beach,  CA,  28-30  Nov  1989 


14.  ABSTRACT 

see  report 


15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 
ABSTRACT 

18.  NUMBER 

OF  PAGES 

19a.  NAME  OF 

RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Same  as 
Report  (SAR) 

7 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


distributions  of  |  +m  >  and  |  —m  >  have  the  potential  to  generate  a  frequency  offset  by  causing 
a  sloping  background  to  the  |  0  >— >|  0  >  transition.  Both  ^-transitions  (Am  —  0)  and  weak  x— 
transitions  (Am  =  ±1)  may  contribute  to  this  background.  Since  the  background  and  the  servo 
gain  slope  have  different  microwave  power  dependencies,  the  resulting  offset  provides  a  means  for 
transforming  power  changes  to  frequency  changes!4! .  Cavity  pulling  can  render  cesium  clocks  sensitive 
to  both  microwave  power  level  and  temperature  changes.  Both  offsets  are  mediated  by  exactly  those 
contributions  which  could  be  safely  neglected  above. 

Apart  from  transient  effects  in  frequency  modulation  schemes,  the  unwanted  signals  can  be  calculated 
by  allowing  the  atoms  to  experience  a  constant  microwave  frequency  during  resonance,  but  one  that 
changes  in  time  before  succeeding  atoms  arrive.  This  preserves  the  conditions  for  standard  microwave 
resonance.  The  transient  conditions,  on  the  other  hand,  require  a  different  description.  This  paper 
describes  an  approach  to  include  the  transient  contributions  to  both  Rabi  and  cavity  pulling.  The 
results  comparing  sine-wave  and  square-wave  phase  modulation  offsets  are  compared  with  experiments 
performed  on  a  commercial  frequency  standard. 

RABI  PULLING 

First  order  time-dependent  perturbation  theory  is  used  to  calculate  the  background  signal  resulting 
from  transitions  which  are  far  off-resonance.  This  approach  can  be  easily  applied  to  any  modulation 
scheme  because  it  allows  the  phase  modulation  to  be  addressed  by  Fourier  decomposition.  Writing 

<  F  =  4,  m  |  M  |  F  =  3,  m1  >=  2 hcn  cos (u>t  +  n<umf  +  0n)  (1) 

n 

for  a  particular  [  m  >->|  m1  >  transition  leads  to  a  general  first  order  transition  probability,  W,  for  a 
single  velocity  given  by 


w  = 


sin(A  +  no>m/2)7"sin(A  +  mwm/ 2)i 


(A  +  rwjm/ 2)  (A  +  rnwm/ 2) 

|cos[(m  -  n)u)m[t  -  T(  1  +  D/L))  +  9m-  6n\  +  cos[(m  -  n)u>m(t  -  TD/L)  +  0m~  0„]}.(2) 


The  notation  follows  [2]  where  2A  represents  how  far  the  microwave  frequency  is  above  resonance,  r 
is  the  transit  time  through  one  arm  of  the  Ramsey  cavity,  T  is  the  transit  time  between  arms  distance 
L  apart,  and  D  is  the  detector  distance  from  the  second  arm.  Velocity  averaging  has  been  presumed 
in  order  to  average  to  zero  the  interference  term  between  interaction  regions  (there  are  no  Ramsey 
oscillations  in  the  wings  of  a  transition!).  The  sin  x/x  terms  arise  from  assuming  perfect  square-wave 
amplitude  variation  of  the  microwave  field  on  passing  through  a  cavity  arm!4!.  The  expression  can  be 
simplified  by  linearizing  these  terms  with  respect  to  frequency  over  the  small  n,m  values  of  interest. 
Then  it  can  be  applied  to  a  particular  modulation  scheme  to  predict  the  fundamental  Fourier  signal 
component  appropriate  to  narrow-band  amplification  and  detection. 


a)  Sine-wave  modulation:  cu(t)  =  u>  +  Ao>coscumt  with  Cn  =  cJn( and  0n  =  0, 


W  _>  ^E-jc  ^2  Arjjcos^m[t  -  T(  1  +  D/L)^  +  cos -  TD/L^  j.  (3) 

This  result  has  been  derived  previously  by  De  Marchi!4!. 
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b)  Square-wave  modulation:  cos[tu£  +  9{t) j ,  with  9(t)  alternating  between  0  and  9  at  frequency 
and  for  which 


cq  =  ccos  9/2,  9q  =  9/2)  Cn  =  (2 c/nx)  sin  0/2,  9n  =  9/ 2  for  n  odd. 


Only  terms  with  n  =  0,  ±1  contribute  and 


W 


sin^'N  d  f  c2  sin2  Ar 


3  A  1  A2 


cosu;m[f  —  T(  1  +  D;  L)  -j-  cos ojm[t 


D/q}. 


(4) 


Aw  times  the  phase  factors  in  Eq.(3)  may  be  interpreted  as  the  instantaneous  frequency  in  the  two  arms 
of  the  Ramsey  cavity.  (2uijn/7r)  times  the  similar  phase  factors  in  Eq.(4)  are  the  Fourier  components 
of  delta  function  responses  from  atoms  residing  in  the  two  interaction  regions  at  the  time  when  the 
step  changes  in  phase  occur.  These  expressions  for  W  have  to  be  averaged  over  velocity  distributions 
and  summed  over  separate  I  m  >— »j  m'  >  transitions  to  complete  the  calculation. 

These  signals  are  nullified  by  a  servo  loop  error  signal,  with  the  resulting  offset  determined  by  the 
servo  gain  slope.  The  gain  slopes  are  determined  for  the  usual  mono-velocity  expressions’13!  which 
must  also  be  separately  averaged  over  the  velocity  distribution  for  the  ]  0  >  — > j  0  >  transition.  All  the 
averaging  procedures  have  been  suppressed  for  the  sake  of  clarity  in  the  following  expressions.  The 
procedures  have  been  discussed  elsewhere^1'. 


(w  -  CUc)a  = 


(Au;)<9  /  c2sin2Ar 
8  A  \  A2 


|  cos u>mT  /  2 cos u>m [t  ■■■-  T(l/2  Dj L)\ 


Tsnr  2ctJi(/3)  coscum[i  —  T(J/2+  D/L)} 


(5) 


(  Au>  \ 

with  j3  =  2( -  j  sinwmT/2. 

V  J 


(w  Wc'jgq  ■ — 


2  lu„,  sin  I 


J3A  ) 


c2  sin2  At 
A2 


cos  wmT /2  cos  |t  —  T(l/2  —  D/L)] 


(6) 


T ( 2 s|n 9 )  sin2  2ct sin comT /2 cos um[i  —  T(  1/2  +  D/L)] 

In  the  limit  of  small  modulation  index,  Aoj  — ►  0  and  (tu  —  ojc ) 3  =  (w  —  ojc)sq.  At  finite  values  for 
Aw,  (cj  —  6u)3  >  (u  —  ojc)Sq.  Exact  values  for  this  inequality  will  depend  upon  the  velocity  averages 
and  the  specific  values  chosen  for  servo  parameters.  However,  the  magnitudes  of  the  offsets  will  remain 
very  similar. 

The  perturbation  expansion  developed  for  W  is  in  powers  of  (c/A).  The  ratio  of  second  order  to  first 
order  terms  is  (c/A)2.  Noting  that  c  rs  2.5  kHz  and  A  rs  20kIIz,  it  is  apparent  that  the  series  is 
rapidly  convergent. 


CAVITY  PULLING 

Cavity  pulling  is  caused  by  the  asymmetry  of  the  microwave  spectrum  about  its  unperturbed  frequency, 
induced  by  a  difference  between  that  frequency  and  the  cavity  resonant  frequency.  The  fact  that 
square-wave  phase  modulation  leaves  the  frequency  unchanged  except  during  the  step  changes  in 
phase  leads  to  a  natural  modeling  approach  in  which  the  atomic  resonance  is  treated  exactly  before 
and  after  the  step  change,  while  the  transition  itself  is  treated  by  perturbation  theory.  Only  those 
atoms  residing  in  the  separate  arms  of  the  Ramsey  cavity  at  the  time  of  the  step  changes  are  affected. 
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If  the  transition  occurs  infinitely  fast,  i.e.  in  a  time  short  compared  to  the  precession  period  of  the 
atoms  in  the  microwave  field,  the  sudden  approximation  can  be  invoked  to  predict  that  the  transient 
itself  will  not  induce  any  atomic  response.  Of  course,  the  step  change  in  phase  will  result  in  a  servo 
signal  from  those  atoms,  which  will  merely  be  added  to  the  usual  servo  signal  arising  from  the  greater 
number  of  atoms  in  flight  between  the  two  arms  of  the  cavity.  However,  the  finite  bandwidth  of  the 
microwave  cavity  prevents  infinitely  fast  changes  in  phase  from  occurring  and  causes  a  breakdown  of 
the  sudden  approximation.  Now  the  transient  frequency  changes  are  finite.  Moreover,  the  positive  and 
negative  frequency  transients  remain  symmetric  only  if  the  microwave  frequency  is  located  at  cavity 
resonance.  It  is  this  asymmetry  which  records  the  cavity  offset  and  induces  a  cesium  beam  signal  at 
the  modulation  frequency.  This  signal  is  derived  from  fast  responses  at  the  beginning  and  end  of  each 
half-cycle  of  the  servo  error  signal  and  may  be  written  as 

—  ^ ^  ^  Ql  sjn2  2cr|(2cr)/tan(2cr)  j  cosc umT /2  cosc jm[t  —  T(l/2  +  D/ L)\  ^ - -  J  .  (7) 

When  this  is  combined  as  before  with  the  servo  gain  slope,  the  offset  (w  —  wc)  can  be  calculated.  For 
a  mono-velocity  beam,  it  reduces  to 

SQl^  co^„r/2| 

w02T  smwmT/2 1/  W  1  Jj\  wo  J  K1 

The  similar  result  for  sine-wave  modulation^  may  be  written 

i^r).  =  11 + Mf>) ~ Jj(/,)}{(2cr)/ta”(2CT)} (^)  (9) 

In  the  limit  Aw  — ►  0,  the  two  expressions  become  equal,  and  cavity  pulling  effects  are  the  same  for 
sine-wave  and  square-wave  modulation.  The  sine-wave  frequency  offset  is  reduced  as  Aw  increases 
from  zero. 


EXPERIMENTAL  RESULTS 

Sine-wave  and  square-wave  modulation  were  compared  by  measuring  the  frequency  of  an  atomic 
standard  for  both  modulation  schemes,  under  as  identical  conditions  as  possible.  This  allows  the 
differences  in  modulation  schemes  to  become  apparent  but  does  not  attempt  to  optimize  performance. 
The  square-wave  modulation  depth  9  was  held  close  to  90  degrees  and  the  modulator  drive  signal 
was  passed  through  a  narrow  band,  unity  gain  filter  and  phase  shifter  to  generate  the  sine-wave 
modulation.  This  produced  a  modulation  index  (Aw/wm)  of  unity  and  allowed  phase  adjustment  to 
maximize  the  signal.  The  modulation  frequency  was  close  to  maximizing  the  gain  slope  of  the  servo 
loop  with  wtoT  »  n  and  /3  «  2. 

Rabi  pulling  was  separately  identified  by  measuring  the  frequency  as  a  function  of  C-field  at  two 
microwave  power  levels;  the  optimum  power  maximized  the  average  beam  current,  and  the  higher 
power  was  3  dB  above  optimum.  Figure  1  shows  the  resulting  fractional  change  in  frequency  with 
power  level  for  both  modulation  schemes.  Cavity  pulling  was  determined  by  measuring  the  frequency 
for  cavity  tuning  on  resonance  and  ±12  MHz  off-resonance.  In  the  data  reduction,  it  was  assumed 
that  the  Rabi  pulling  was  independent  of  cavity  tuning,  and  the  cavity  pulling  was  independent  of 
C-field.  Three  cavity  settings,  two  power  levels  and  two  modulation  schemes  lead  to  a  total  of  twelve 
curves  of  frequency  versus  C-field.  The  Rabi  pulling  could  be  reduced  to  a  single  curve  at  optimum 
power  by  the  following  relationships: 
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(a>  -  u)c)3q 
and  (a;  —  wc)(P  +  3dB 


(10) 

(11) 


=  0.94(a>  —  u>c)s  for  both  power  levels, 

=  1.85(cu  —  uc)(P)  for  both  modulation  schemes. 

This  curve  is  shown  in  Fig. 2  and  illustrates  that  Rabi  pulling  effects  are  small  in  this  standard.  Zero 
frequency  offset  may  be  arbitrarily  chosen  for  the  highest  C-field,  but  is  trucly  zero  only  at  even  higher 
fields.  The  ratio  0.94  between  square-wave  and  sine-wave  amplitudes  is  closer  to  1  than  predicted 
from  Eqs.(5)  and  (6).  In  the  absence  of  velocity  averaging,  this  ratio  is  predicted  to  be  di(2)  =  0.58. 
Also,  the  ratio  between  high  power  and  optimum  power  is  predicted  to  be  2.5,  slightly  higher  than 
the  measured  value  of  1.85.  It  is  reasonable  to  anticipate  that  velocity  averaging  will  improve  the 
agreement  between  theory  and  experiment.  The  approximately  linear  dependence  on  power  arises 
because  the  off-resonance  transition  probabilities  increase  linearly  with  power,  whereas  the  resonant 
transitions  are  already  optimized. 

Table  I  lists  the  measured  offset  frequencies  versus  cavity  tuning.  Systematic  errors  of  unknown  origin 
exceed  the  standard  deviation  (2  X  10-13)  in  the  frequency  measurements.  Zero  offset  was  arbitrarily 
chosen  to  be  that  for  the  sine-wave  modulation  case  at  optimum  power  and  cavity  tuned  on  resonance. 
Eqs.(8)  and  (9)  would  predict  a  mono-velocity  offset  ratio 

(w  —  (ue)aq  2Ji(2)  _  , 

(cu-Wc)s  1  +  J0(2)-  J2(2)  ‘  1  J 

The  measured  ratio,  obtained  by  taking  the  difference  of  the  offsets  at  ±  12MIIz,  is  1.3  at  both  power 
levels.  The  power  dependence  of  the  offsets  is  contained  in  the  usual  expression  (2cr)/tan(2cr)  which 
is  zero  for  4cr  =  7r  and  rises  to  —1.7  as  the  power  is  raised  3  dB.  Clearly,  the  optimum  power  did  not 
produce  zero  cavity  pulling,  but  the  increase  by  a  factor  2.7  on  raising  the  power  3  dB  would  imply 
that  the  nominal  power  was  within  1.6  dB  of  that  necessary  for  zero  offset. 


CONCLUSIONS 

Our  measurements  show  that  sine-wave  and  square-wave  phase  modulation  yield  similar  values  for 
Rabi  and  cavity  pulling  effects  in  cesium  frequency  standards.  The  result  is  consistent  with  the  theory 
presented  in  this  paper.  It  predicts  that,  in  the  limit  of  small  frequency  modulation  depth,  sine-wave 
modulation  produces  the  same  frequency  offsets  as  square-wave  modulation.  At  higher  modulation 
depth,  however,  Rabi  pulling  tends  to  be  larger  and  cavity  pulling  smaller  for  the  sine-wave  modulation 
scheme.  Exact  predictions  will  require  knowledge  of  specific  servo  parameters  and  averaging  over  beam 
velocities. 
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TABLE  I. 

Fractional  frequency  offsets  (parts  in  1012)  determined  from  measurements 
at  two  microwave  power  levels  for  both  sine-wave  and  square-wave  modula¬ 
tion.  The  cavity  was  tuned  to  the  cesium  frequency  Fq  and  to  Fort  12MHz. 
The  cavity  pulling  is  normalized  to  the  frequency  for  sine-wave  modu¬ 
lation  at  optimum  power  (maximum  dc  signal)  and  cavity  tuned  to  Fq. 
Standard  deviation  =  3  X  10" 13. 


F0  +  12  MHz 

Fq 

Fq  -  12  MHz 

optimum  power 

sine-wave 

-0.9 

0 

3.6 

square-wave 

-3.6 

-1.2 

2.2 

optimum  +  3dB 

sine-wave 

-3.6 

-0.2 

8.5 

square-wave 

-8.2 

-1.0 

6.9 
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Figure  1.  C-field  dependence  of  the  fractional  change  in 
frequency  caused  by  a  3  dB  increase  in  microwave  power 
above  optimum  (maximum  dc  signal)  for  sine-wave  and  square 
wave  phase  modulation. 


30  40  50  60 


7FFMAN  FRfcQUtNCV  (KH/) 

Figure  2.  C-field  dependence  of  Rabi  pulling  measured  on 
a  cesium  frequency  standard  with  sine-wave  modulation. 
Rabi  pulling  is  assumed  to  be  those  changes  in  frequency 
which  depend  upon  C-field. 
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